Perfect Impedance-Matched Isolators and Unidirectional Absorbers 
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A broad-band reflectionless channel which supports unidirectional wave propagation originating 
from the interplay between gyrotropic elements and symmetrically placed gain and loss constituents 
is proposed. Interchange of the active elements together with a gyrotropic inversion turns the 
same structure to a unidirectional absorber where incoming waves from a specific direction are 
annihilated. When disorder is introduced asymmetric Anderson localization is found. Realizations 
of such multi-functional architectures in the frame of electronic and photonic circuitry are discussed. 
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Understanding and controlling the direction of wave 
propagation is at the heart of many fundamental prob- 
lems of physics while it is of great relevance to engineer- 
ing. In the latter case the challenge is to design on- 
chip integrated devices that control energy flow in differ- 
ent spatial directions. Along these lines, the creation of 
novel classes of integrated photonic, electronic, acoustic 
or thermal diodes is of great interest. Such unidirectional 
elements constitute the basic building blocks for a vari- 
ety of transport-based devices like rectifiers, circulators, 
pumps, coherent perfect absorbers, molecular switches 
and transistors [IHI]. 

The idea was originally implemented in the electronics 
framework, with the construction of electrical diodes that 
were able to rectify the current flux pQ. This significant 
revolution motivated researchers to investigate the pos- 
sibility of implementing the notion of "diode action" to 
other areas of physics ranging from thermal [5] and acous- 
tic transport [B] to optics [2]. Specifically in optics, uni- 
directional elements rely almost exclusively on magneto- 
optical (Faraday) effects. However, at optical frequencies 
all nonreciprocal effects are very weak resulting in pro- 
hibitively large size of most nonreciprocal devices. Al- 
ternative proposals include the creation of optical diodes 
based on nonlinearities [H [7H1] ■ These schemes how- 
ever suffer from serious drawbacks and limitations since 
the rectification depends on the level of incident power 
or/and whether the second harmonic of the fundamental 
wave is transmitted or not. Obviously, in this later case 
the outgoing signal does not have the same characteristics 
as the incident one. Moreover non-linearities can result 
in enhanced reflection, which significantly compromises 
the diode performance. 

In this Letter, we propose a physical setting (see Fig. 
[I]) which acts as a perfect one-way valve i.e. a chan- 
nel along which waves propagate in only a single direc- 
tion, with zero reflection, a perfect impedance match- 
ing isolator. This is achieved by employing an inter- 
face between a gyrotropic element and two active con- 
stituents, one with gain and another one with a bal- 
anced amount of loss. In contrast to standard parity-time 



('PT)-symmetric systems, first proposed in the frame- 
work of optics by Christodoulides and colleagues [TO] , our 
structure shows a new type of generalized 'PT-symmctry 
which allows for non-reciprocal transport. Furthermore, 
an exchange of the active elements together with an in- 
version of the gyrotropy leads to unidirectional absorp- 
tion where an incoming wave entering the structure from 
one side is completely annihilated. These features can 
be observed over a broad range of frequencies when ID 
chains of such elements are considered. Finally, the pres- 
ence of imperfections results in destructive interferences 
that are sensitive to the direction of propagation, a phe- 
nomenon refer to as asymmetric Anderson localization. 
Below we present realizations of such multi-functional ar- 
chitectures, both in the frame of photonics and electron- 
ics circuitry (EC). The later framework has been proven 
recently [TTJ extremely useful for the investigation of the 
transport properties of 'PT-symmetric systems. 

The photonic structures that we consider consist of 
a central magneto-optics layer sandwiched between two 
equally balanced gain and loss birefringent layers which 
allow for coupling between the x — y polarizations (see 
Fig. [TJi). The role of the magnetic element is to in- 
duce magnetic non-reciprocity which is associated with 
the breaking of time reversal symmetry. Breaking 
time-reversal symmetry is not sufficient to obtain non- 
reciprocal transport - the theory of magnetic groups 
shows that the absence of space inversion symmetries is 
also required. This is achieved with the use of the two 
active birefringent layers [12] and an anisotropic Bragg 
mirror which has missalignment with the active layers; 
therefore it does not allow for coupling between the two 
polarization channels. 

The equivalent EC (Fig. [T]b) comprised of two pairs of 
mutually inductive, M = /j,L, coupled LC oscillators, one 
with amplification (left column of Fig. [lb) and the other 
with equivalent attenuation (right column of Fig. [lb). 
The loss imposed on the right half of the structure is a 
standard resistor R. Gain imposed on the left half of the 
EC, symbolized by -72, is implemented with a negative 
impedance converter (NIC). The NIC gain is trimmed to 
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oppositely match the value of R used on the loss side, set- 
ting the gain and loss parameter 7 = The un- 
coupled frequency of each resonator is cj = 1/VLC. The 
pairs are coupled with each other via a gyrator [TB . This 
element immitates the role of the magnetic layer used in 
the photonic set-up. The gyrator is a lossless two-port 
network component with an antisymmetric impedance 
matrix Z g = —Zj connecting the input and output volt- 
ages V = (V n ,V m ) T and currents I = (I n ,Im) T associ- 
ated with ports n and m as 
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where R g = j3 1 y ^ and /? is a dimensionless conduc- 
tance. Eq. ([!]) is invariant under a generalized time 
reversal operator T which performs a combined time- 
inversion (t — > —t) together with a transposition of Z g . 
The anti-linear operator T is equivalent to the transfor- 
mation t —> —t; (3 —> —p. Despite the fact that a gyra- 
tor breaks the time-reversal symmetry, it does not alone 
leads to non-reciprocal transport (in a similar manner 
that a magnetic layer is not sufficient to create asym- 
metric transport in the optics framework). 

Finally, the two (upper/lower) propagation channels 
that are supported by the EC structure of Fig. [TJd imitate 
the x/y polarization channels of the photonic set-up. In 
the latter case, the BG act as a polarization filter for inci- 
dent waves with frequencies chosen from the pseudo-gaps. 
These pseudo-gaps are a consequence of the anisotropic 
refraction index along the x and y polarization channels. 
The analogous effect can be achieved in the EC case by 
coupling only the upper (or lower) channel to a transmis- 
sion line (TL) with impedance Z . 

Although below we exploit the simplicity of the elec- 
tronic (lump) circuitry framework in order to guide the 
physical intuition and derive theoretical expressions for 
the transport characteristics of these structures, we will 
always validate our results in the realm of photonic cir- 
cuitry via detail simulations. 

At any point along a TL, the current and voltage deter- 
mine the amplitudes of the right and left traveling wave 
components. The forward V^, R and backward V/~ m wave 
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amplitudes, and Vl/r and Il/r the voltage and current 
at the left (L) or right (R) TL-EC contacts satisfy the 
continuity relations 
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Note that with this convention, a positive lead cur- 
rent flows out of the right circuit, but into the left cir- 
cuit, and that the reflection amplitudes for left or right 
incident waves are defined as t^r = V L f R /V L t R and 

Il,r = V^ l / V l,r respectively. 

Application of the first and second Kirchoff 's laws at 
the TL-EC contacts allow us to find the current /voltage 
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FIG. 1: (Color online) A "PT-symmetric photonic struc- 
ture. Red and green layers indicate gain and loss. The grey 
layer indicates a magneto-optical material. The orange ar- 
row indicates the direction of the magnetic field (gyrotropy) . 
Light/dark blue layers indicate a Bragg grating which is 
missaligned with the active layers, (a) A perfect impedance- 
matched isolator where reflection is vanished, (b) The same 
structure can also act as a unidirectional absorber if we ex- 
change the gain/loss layers and invert the gyrotropy. (c) The 
equivalent electronic circuit. R g indicate the gyrator. 



wave amplitudes at the left and right contact. We get 



AlV L - P(Vi - V 2 )] - ujV l + 
ihVR-i3(V 1 -V 2 )}+uV R - 
i[-yV 1 + P(V L -V R )]-uV 1 + 
i[-rV 2 + /3(V L -V R )}+uV2- 
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where rh = 1 — /i 2 , Vj.,2 is the voltage at the lower 
channels (see Fig. [iji), 7 is the gain/loss parameter, and 
r] = \j 'L/C /Zo is the dimensionless TL impedance. Here, 
the dimensionless wave frequency uj is in units of loq. 
Equations ([3| are invariant under a combined parity V 
(i.e. L o R, 1 O 2) and f (i.e. i -B- -i,/3 -H- -(3) 
transformation (we recall that the time-reversal t — > — t 
is equivalent to conjugation). 

For the EC structure of Fig. [lb, the transmission and 
reflection amplitudes in terms of the gyrotropic and dis- 
sipative conductances j3 and 7 are derived using Eqs. ^ 
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where the functions /(...), g(-) and h(...) are given 
at the supplement and satisfy the symmetry relations 

/(w;/3,7) = /(w;-j0,-7)j A 7) = g(u;-(3,j) and 
ft.(w;/3,7) = h(ui; \/3\, Consequently, these symmetry 
relations can be translated to the following relations for 



the transmission and reflection amplitudes for any fre- 
quency u 

t a (P,'y)= t R (-p,-i) =t L (-p,i) = t L (p,-j) 
r R (/3,'y)= r R (-P,<y) = r L (P, -7) = r L (~(3,-j). 

(5) 

The corresponding left (L) and right (R) transmitancc 
and reflectance are then defined as Tl/r = \t L /R\ 2 an d 
Rl/r = Vl/r? respectively. 

The ratio between the left and right transmission takes 
the simple form 



t_L_ 
tR 



= 1 



2 7M 



p - - /3(1 - M 2 )^ 



2 ' 



(6) 



which establishes the nonreciprocal nature of Vl~- 
symmetric transport. Note that 7 = results in recip- 
rocal transport irrespective of the value of p. The same 
conclusion is reached for P = 0, corresponding to stan- 
dard 'PT-symmetric structures. Therefore we conclude 
that non-reciprocal transport is achieved due to the the 
interplay between gain/loss elements (i.e. 7) together 
with gyrotropic constituent p. 

Furthermore one can show that the left (right) trans- 
mission tuR\ and reflection tur) amplitudes satisfy the 
relations 



Re[ t L r* R } = 0, Re[ t R r* L ] = 0, 
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which again implies that in general transport is non- 
reciprocal. In the photonic case the transmission and 
reflection amplitutes becomes 2x2 matrices. In fact, 
Eqs. ([7]) are consequences of generalized unitarity rela- 
tions satisfied by the scattering matrix S [HI [T3] . 

In Fig. [2^, we report representative cases of the re- 
flected and transmitted signals for a left and right inci- 
dent waves. A striking feature of these plots is the fact 
that at specific frequencies uj p the transmittance from one 
side can become zero, while from the other side is high. 
At these specific values of (w p ,7 P ,^ p ), the EC acts as a 
perfect isolator. For example, the requirement for zero 
transmittance for a left incoming wave t L = leads to 
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while the frequency where tR = is simply w p (— 7p , P p ). 
At these frequencies Eq. ^ collapses to the relation 
RlRr — 1, he. the reflection from one side is enchanced 
while from the other side is suppressed. 

Depending on the reflectance properties at u) p , we now 
distinguish two different types of perfect isolation: (a) In 
the first case, the reflection of an incident wave entering 
the sample in the direction that we have non-zero trans- 
mitivity, is zero. For example, the transport characteris- 
tics of a right incident wave would be T R ^ and Rr = 
while, the same wave entering the circuit from the oppo- 
site (left) direction results in Tl = and R L 7^ 0. We 
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FIG. 2: (Color online) (a) Transmissions and reflections 
Tl/r', Rl/r of one EC unit for a left/right incident wave. The 
red circle indicates a UA while the blue triangle indicates an 
IMI frequency for a right incident wave; (b) A broad-band UA 
for a right incident wave for a periodic ID array of N — 10 EC 
units. We have checked (not shown) that changing /3 — > — (3 
we get an IMI for the same frequency window (see text); (c) A 
broad-band UA for an a;-polarized right incident wave in the 
case of a periodic photonic structure; (d) Changing /?—>—/? 
leads to an IMI behavior (see text) for the set-up of subfigure 
(c); (e) Additional change of 7 — > —7 leads to UA for a left 
incident wave. 



refer to this case as impedance-matched isolator (IMI) as 
the transmitted signal does not experience any reflection 
at the TL-circuit interface, (b) In the second case the 
signal is completely absorbed (i.e. neither transmitted 
nor reflected) from one direction, while it is transmitted 
(and partially reflected) from the other. We refer to this 
case as unidirectional absorber (UA) . An example of such 
modes associated with the former case is marked with a 
(red) circle in Fig. [2^, while a representative mode for 
the latter case is indicated with a (blue) triangle in Fig. 
[2^,. Based on the symmetry relations of Eq. ^ we de- 
duce that a change from a UA to a IMI behavior occurs 
by simply performing the transformation P — > — p. In- 
stead a simultaneous exchange of gain/loss elements i.e. 
7 — > —7, will result in the same behavior (UA or IMI) 
but for an incident wave entering the structure from the 
opposite direction. 

In Fig. [2)3 we report the left /right reflectances and 
transmit ances, for a ID network consisting of Af iden- 
tical EC units. We observe that a perfect isolation is 
now extended over a broad-band frequency window (see 
highlighted area). A simple theoretical understanding for 
the isolation action can be achieved by realizing that in 
this case t% /t% = det 1% = (detM)^ where M is the 
transfer matrix of each EC unit and Ma/ is the transfer 
matrix associated with the total array. A Taylor expan- 
sion of the right hand side around a frequency lo p of an 
individual unidirectional unit requires knowledge of the 



-th order derivatives at w D i.e. 



d n del. 



for 



n = 1, • • ■ ,M — 1. These vanishing derivatives apprecia- 
bly flatten the transmitance at a neighborhood around 
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ujp, effectively creating perfect isolation for a broad-band 
frequency domain, in agreement with the simulations 
shown in Fig. [2p. Moreover an IMI behavior is evident 
from our data. We have also checked that the previ- 
ous expectations associated with the /3, 7-symmetries are 
confirmed in the case of ID chains as well. The same be- 
havior is also present for the equivalent periodic photonic 
crystal consisting of Af basic unidirectional units as can 
be seen from Figs. [2j>e. 

Finally, we investigate the effect of imperfections along 
ID chains with local PT-symmetry. We assume that in 
the case of ECs the disorder is introduced on the capac- 
itances C which are independently distributed around a 
typical capacitance Co with a uniform probability dis- 
tribution in the interval [CijCy. In the absence of 
any gain/loss elements, such ID chains exhibit the phe- 
nomenon of Anderson localization. This phenomenon 
was predicted fifty years ago in the framework of quan- 
tum (electronic) waves by Anderson |16j and its existence 
has been confirmed in recent years by experiments with 
classical [TTH2"U] and matter waves [5TJ [55] . 

Formally, the Anderson localization can be quantified 
via the so-called localization length £ which is defined as 
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where (• • • ) denotes an ensemble average over disorder 
realizations. In reciprocal disordered systems, the local- 
ization length satisfies the relations £l — £,r- We find 
instead that our set-up allows for asymmetric localiza- 
tion with £l ^ £r. This is demonstrated in the upper 
inset of Fig. [3] where we report the scaling behavior of 
(lnXi^) for an ID chain consisting of Af disordered EC 
units, each respecting locally a PT-symmetry. Straight- 
forward algebra leads to the following expression for the 
difference between the two localization lengths 
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where A^ 1 = a±(/s) = /3(1- K (1-m 2 V 2 )±7M 

and K12 — C12/C0. In the main panel of Fig. [3]we com- 
pare the theoretical predictions of Eq. ( 10 1 together with 



the numerical results from the transfer matrix simula- 
tions of the disordered chain. 

The asymmetric localization £j, 7^ can be also ob- 
served in the equivalent photonic crystal consisting of Af 
random PT-symmetric units similar to the one shown in 
Fig. [TJi. The disorder was introduced in the real part of 
the permitivity of the active elements in such a way that 
in all cases the individual units were preserving a local 
PT-symmetry. The numerical evaluation of (ln(2Tft) p ) e 
is performed using transfer matrices and the results are 
shown in the lower inset of Fig. [3j In this case, an addi- 
tional average over different polarizations was performed 




FIG. 3: (Color online) Upper inset: The scaling of {T^rm} 
versus the system size N for a disordered ID chain with ran- 
dom PT-symmetric EC units. Main panel: The difference 
A£ _1 between left and right localization lengths for various 
Jt2 («i is kept fixed). The numerical results using transfer 
matrices (blue squares) are compared with the theoretical ex- 
pression (solid line) of Eq. Lower inset: The scaling of 
(Tlir)) versus the system size M for the equivalent disordered 
photonic hetrostructure consisting of PT-symmetric units. 



in order to guarantee that the phenomenon is polariza- 
tion independent. 

In conclusion, we have demonstrated both theoreti- 
cally and numerically the dual behavior of PT-symmetric 
(photonic and electronic) circuits as perfect-impedance 
matched isolators or unidirectional absorbers. One- 
dimensional arrays consisting of such PT-symmetric 
units can demonstrate these dual properties for a broad- 
band frequency window. Finally, we have shown that dis- 
ordered ID arrays with local PT-symmetry, show asym- 
metric (left /right) interference effects that can lead to 
asymmetric Anderson localization. It will be interest- 
ing to exploit the simplicity of ECs in order to inves- 
tigate experimentally such phenomena and extend the 
analysis towards directional chaos and localization where 
a localization-delocalization phase transition was pre- 
dicted even for ID arrays [23"H2"5] , We expect that fab- 
ricating such multi-functional components in the frame 
of electronic complementary metal oxide semiconductor 
(CMOS) circuits [Mj and the frame of Indium Phosphide 
Photonic (InP) monolithic Integrated Circuits [57] will 
occur the years to come and the findings of the optics 
community will pollinate the silicon photonics and InP 
communities. 
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Supplementary Material for Transmis- 
sions/Reflections AMPLITUDES [EQ. (4)] 

The transmissions and reflections amplitudes can 
be expressed using the following three functions /(7, /?), 
#(7, /?), and h(j,(3) which are defined for any frequency 
lo as 

, f(-r,P) . . /(~7,fl . r 9(1,(3). 

tR= — 7 — ' *i = r ? r L — — rr — , 

h n in 

where / and g are real functions such that 

/(7,/3) = /(-7,-/3) 9{1,P)=9(1,-P) 

To express these three functions the following sho: 
notations are used: 

M^(l- M 2 ) 
K = M ( 7 2 + uj 2 ) - 2 

A = 2M (M ( 7 2 - 2^ 2 ) - 2) 
B = M (8/3 2 - 4 7 2 + 7 4 M + 2) + 4 
C = 2 ( 7 2 (V + 1) - 2^ 2 M - 2) 
P = M 2 w 8 + Ad 6 + Blu a + Clu 2 + 1 



Using the above expressions we get after straightforward 
algebra, 



in 



rthand 



/( 7) /?) = 4/3?7Ma; 3 (^Mw 2 - /3 + 7/1) 
5 ( 7 ,/3)=£ + P + £ 
h = R + i(P-D) 



D = tfu? (Mlu 2 K + 1) 

E = 2777a; 2 ( M 2 + Muj 2 K + 1) 

i? = -2 V oj (Moj 2 - 1) (lj 2 (-2M/3 2 + K) + l) 



